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MINITHACK POSTDETECTION 


BANDWIDTH CONSIDERATIONS 

T. J. Grenchik 
C. W. Murray, Jr. 

ABSTRACT 

A description of the Minitrack system and definition of the measurement 
made by the system is presented. 

Design considerations relating to the phaselock loop response of the pro- 
posed postdetection filter are discussed. Under an assumed set of "worst case" 
orbital dynamics for the satellite (100 statute mile height circular orbit with a 
useful 10® overhead arc), a second order phaselock loop with a 10 Hz bandwidth 
is recommended for the filter. 

The "gain" due to least squares smoothing of the Minitrack data at the out- 
put of the postdetection filter is calculated for a number of filter designs. For- 
the 10 Hz equivalent double pole filter above, a sampling rate from 10 to 20 
samples per second is recommended. It is shown that by increasing the rate 
above 20 samples per- second the "gain" due to smoothing increases by approxi- 
mately a db (decibel), while aggravating the data handling problem. 
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SUMMARY 


The Minitrack system and its measurement are described and several po- 
tential error sources related to the system measurement are presented. The 
major source of error listed is the uncertainty in delay (fixed or time varying) 
in each of the antennas of the antenna pair prior to signal combining (or equiva- 
lently signal cross- correlation). From this major error source, a minimum of 
10 arc seconds pointing error is postulated as the Minitrack system resolution. 

A "worst case" orbit is assumed in order to establish rates of change of 
pointing angle which define the postdetection filter parameters. From a circu- 
lar orbit of 100 statute miles height above the Minitrack station, there follows 
the requirement for a 10 Hz bandwidth, 2nd order phaselock loop filter design. 
Other parameters such as signal detection method and antenna gain remain fixed 
in this analysis. Signal to noise considerations then require adjustment of satel- 
lite effective radiated power. 

The effect of least squares smoothing of the data at the output of the post- 
detection filter is investigated for a number of filter designs and bandwidths. 
Again a "worst case" orbit is assumed with the satellite passing directly overhead 
and the total arc of visibility taken as 10® (±5® from zenith). The "gain" in sig- 
nal to noise ratio due to smoothing is defined and is seen to be a function of the 
satellite dynamics, the total time of visibility, the filter design and bandwidth. 
It is shown that the choice of a filter design cannot be based on the "gain" achieved 
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by smoothing. Once the design is fixed the least squares analysis aids in the 
selection of a sampling rate. 

For the 10 Hz equivalent double pole filter recommended above,' a sampling 
rate of 10 to 20 samples per second is recommended. This particular range of 
sampling rates will provide a sufficiently high ’'gain” while keeping the data 
handling at a minimum leyel. 
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MINITRACK POSTDETECTION 
BANDWIDTH CONSIDERATpNS 


INTRODUCTION 

The Minitrack system measures angles (direction cosines) to ah orbiting’ 
spacecraft with reference to a ground station coordinate system by an interferom- 
eter technique. The system has been successfully used for tracking near-Earth 
satellites (reference 1) where angle rates are appreciable as compared to satel- 
lites at lunar distances and beyond. 

This paper considers system design parameters' for the postdetection filter 
of the Minitrack system and' discusses the "gain” due to least squares smoothing 
of data at the output of the filter. The orbital djmamics and geometry which -are 
assumed can be seen in Figure 1. 

The satellite located at point S is in a" circular orbit at an altitude h ^ above 
the earth and. is moving in a clockwise direction. The line of sight from the 
tracking, station at point T to the satellite makes an angle Q (interferometer 
angle) with the local tangent. The central angle is a and the angle between line 
ST and SO is /3 . The system measurement is proportional to D cos d where D 
is the length of the interferometer baseline. 
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From the changing geometry of Figure 1, design requirements are imposed 
upon the postdetection filter which define the filter shape and bandwidth. Selection 
of Minitrack system data rate follows this filter definition and it is shown that an 
optimum range of sampling rates exists that is practical with respect to 
Goddard/Minitrack Station- data handling capabilities. 

SYSTEM DESCRIPTION 

Figure 2 shows a simplified block diagram of the Minitrack system. In this 
diagram a pair of antennas, separated by a baseline length, D, receives a signal 
sini// (t), where xp{t) is the time variant phase function of the incoming signal. 

The knowledgeable reader will immediately note that the ambiguity resolution 
system has not been included. This was done to simplify the description. For 
the reader without a Minitrack system background, References 2,3, and 4 are 
recommended. The pair of fine antennas imparts a phase change U [<9 (t)] and 
V [^(t) ] to the signal observed at the respective antennas. Because of the 
angle the line-of-sight angle with respect to the local horizontal, a point on the 
wavefront, sin i/;(t), arrives at FINE ANTENNA ”B" a time At later than at FINE 
ANTENNA ”A". The difference in path length, 0, is related to the line-of-sight 
angle, & , by:^ 

$ = D cos 0 

^ Refraction has. not been considered in this paper. See reference 12. 
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and the time of arrival difference is given by: 

- T, = At (2) 

with the time of arrival at antenna "B”, T the time of arrival at antenna "A", 
$ and D in meters, Q in radians. At in seconds, and c, the speed of light, in 
units of meters/second. The measurement made by the system is related to 
the time of arrival difference At. How this occurs is seen by following the antenna 
inputs through the simplified block diagram in Figure 2. 

At antenna "A", at some instant of time the incoming signal is sin 
At that same instant of time, the observed signal at antenna "B” is: 

D cos (3) 

For example, if the input at antenna *’A” is a constant frequency sinusoid, then • 

sin ip (t) = sin (2tt f ^ t) (4) 

where f ^ is the frequency of emission from the satellite in Hz« This assumption 
neglects Doppler effects which are included in later sectionSo At Antenna "B"; 
sin (t - At) = sin [2 t 7 f ^ (t - At)] 

/ - 2 7T f . D \ 

sin ip (t - At) = sin f 2 t cos Oj \°) 


sin ip (t - At) = sin 


K 
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Examine the term 27 rf^ D cos 0 /c. The quantity f^c is equal to l/k^ where 
is the wavelength of the transmitted frequency and is assumed a constant. 


27Tf^Dcos6> 27 t_ . 

— p cos 6 

c k 


(7) 


The baseline separation, D, is approximately equal to 50\^at a frequency f^ of 
136 MHz, the nominal operating frequency of the Minitrack system. Assume for 
the moment that D is exactly equal to 50 Then 


D cos 0 = 2 7T (50) cos 6 (8) 

Since cosine $ can vary between +1 and -1 for 0° < ^ < 180°, the effective phase 

difference {2n/k )D cos d between the two antennas "A" and "B" can vary be- 
1 

tween +50 (2 tt) and -50 {2tt) radians. The Minitrack system (and any interferom- 
eter system operating with a narrow spectrum emission) cannot resolve the 
multiples of 2 7 t with only the two antennas, fine antenna "A” and fine antenna 
"B”, spaced apart many wavelengths of It should be noted that other pairs of 
Minitrack antennas, spaced smaller distances apart, resolve the integral number 
of 2 7 T radians change, and the pair of fine antennas measures the excess angle 
remaining. 

For the moment assume no restriction on the phase fonction (t). At 
antenna "A” output we have 
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sin {ip (t) + U (t)]'} 


( 9 ) 


where U [^(t)] is the phase change to i/f(t) caused by the phase pattern of 
antenna "A" looking in the direction 0 , Of course, 6 is a time varying function. 
Similarly at antenna "B” output we have: 

sin {0 (t - At) + V [0 (t - At)]*} (19) 

but note here the phase change at antenna "B", V[^] differs from the phase 
change U [i9] because of the time of arrival difference. At. 

A reference frequency of 100 Hz, f ^ , is mixed with the antenna signal 
to produce a signal at point B in Figure 2: 

sin {xjj (t - At) +V [0 (t - At)] + 2 7T t} (H) 

No change is made to the antenna *’A" signal and at point A in Figure 2, we have, 

sin {./.(t) +U [0 (t)] } 

The signals at points B and A are summed: 

sin {0 (t - At) + V (t - At)] -t- 2 77 t} + sin {xfj (t) + U [0 (t)] } = 

2 cos i {'2 77 t + </; (t - At) - 0 (t) + V [0 (t - At)] -U [0 (f)] } 

xsini{2?7 t+i//(t-At) + t//(t) + V [^(t-At)] +U [(9(t)]} (13) 
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The appearance of this signal is shown in Figure 3 as "Summer Output." From 
Equation (13), the envelope of the summed signals B + A is seen to carry angle of 
arrival information. In Figure 3, At was made zero and V [0] - = U [^] . The 
envelope is amplitude detected in a half wave rectifier to recover (see Figure 3) : 

cosi{ 27 T t -At) +V [e(t-At)]-U [i9(t)]} 

The postdetection filter output (the subject of this paper) recovers 

cos {277 t - At) -t/;(t) +V [0 (t - At)]-U [0 (t)]-} (15) 

In Equation (15) the nominal frequency of the filter output is 100 Hz, the frequency 
of the reference signal; the if; {t -At)- ifj{t) phase term carries (ambiguous) in- 
formation about the angle of arrival 0 ; and the V [0{t -At)] - U [0 (t)] represents 
the phase term- resulting from the antenna phase perturbations. In theory, these 
antenna phase perturbations are predetermined by aircraft flyover calibrations. 

The significance of Equation (15) is the conversion of the phase difference 
between the two antennas i// (t - A t) - t/f(t), from a phase difference at a nominal 
frequency of 136 MHz, to a phase difference at a frequency of 100 Hz. This 
amounts to a scaling of time of approximately 136 MHz/lOO Hz, or exactly f^ /f^ 
for a constant frequency f ^ . 
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At the postdetection filter output, for each negative going zero crossing of 
the filter ed/detected signal, a stop piilse is generated to control the. stop time of 
a counter which counts a frequency of 100 kilohertz. The start pulse is generatea 
by eactf negative going zero crossing of the reference frequency. Hence the 
Minitrack measurement, , is the number of cycles of f^ , 100 kHz, which 
occur between the negative going zero crossing of the 100 Hz reference frequency, 
f^ , and the succeeding negative going zero crossing of the detected signal (in the 
time interval, S seconds) . See Figure 4 for a pictorial representation of this 
measurement and its ambiguity. 

In Equation (15) , define <f) (t) as the following quantity: 

<P(t) = - At) + V (t - At)] -U [e (t)] } (16) 

Now assume that:^ 

U [e (t)] = V (t - At)]- (1'^) 

The effect of this assumption will be discussed under filter response to satellite 
dynamics. Then 

^(t) = {./.(t - At) -./i(t)} (18) 


^ In present practice, the sample rate is controlled by deleting pairs of start/ stop pulses. For example 
a 5/ sec sample rate would utilize each 20th pair of start/ stop pulse to control the counter. For all 
other pairs the counter would be inactive, 

^ This is equivalent to saying that V [ ^ (t — At) ] -U[0(t)]=K, {&) where K {&) is an exactly 
known precalibrated function. 
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When Figure 4 is examined, it becomes apparent that^(t), the total phase 
difference, is not utilized for the system measurement. Instead, the remainder 
of (t) after division by 2 tt radians, determines the measured quantity. Define 
a quantity 


modulo 2tt (19). 

which is the division of ’^(ty by 27 r radians and the remainder is the result, 

S(tY, .’ Hence the-input to the stop- pulse’ generator is, 

'Mod 2'JT 


cos [2 7T e, ^]- 

and negative going zero crossings (stop pulses) occur at 




k = 0, 1. 2,,. 


(4k + l)|-<^(t)„,, 2. 

2 77 f ^ 


( 20 ) 


( 21 ) 


( 22 ) 


where t - are the negative going zero crossings of the signal. 


The input to the start pulse generator is cos (2 t 7 f ^ t)and start pulses occur 


at: 


k = 0, 1, 2 
(4 k + 1) ^ 

t '= 

. 2 77 . 


(23) 


(24) 


where t are the negative going zero crossings of the reference. 


8- 



The difference in the time of negative going zero crossings, is Equation (22) minus 
Equation (24), defined by the symbol S in seconds ( S 5 ^ At, see Figure 4 for illustra- 
tion, when f is a constant). 


2 7T 


(25) 


S = - 




2 Tr f 


[<P(t) -0(t -At)]^ 

2 7T f 


(26) 


In reference to Figures 2 and 4, it can be seen that 8 in seconds is: 


S = 


(27) 


where C ^ is the number of cycles, or counts, of the counting frequency, f^ , and 
f^ = 100 kHz. 

c 


Assume that the input to antenna "A" is a constant frequency sinusoid, that 
is Equation (4): 

sin v/j (t) = sin (2 77 f ^ t.) 
and the input to antenna "B" becomes: 


sin (t - At) 


/' 2 77 f D cos 0 \ 

3ln(^27r t j 


* 9 


( 6 ) 



Further assume Equation (1*7) and see footnote -3;- 


U 10 Ct)] =Y [0 (t - At)] 
With these assumptions. Equation (26) becomes: 


( 17 ) 


a3 


[2 w t - ^2 7T f, t - 2 ir 


2 7T 


(28) 


Cn 1 T2 -n" f t 1 

; 1.D cos d\ 

f - 2 7T f L c X 


(29) 


From Figure 4, it may be seen that for a constant frequency sinusoid, 


— At = K^+ 5 


(30) 


where K is the ambiguily number determined by the closer spaced Minitrack 
antenna pairs. From Equation (2) 


cos ^ = - A t 
D 


(31) 


^ C TK ^o"| 


(32) 
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For later reference, it is appropriate here to examine magnitudes of typical 
error sources in the measurement. The first source of error to be shown is 
the quantization error of the measurement counter. From Equation (32), 


■ sin d AQ ~ — — 




(33) 


Baseline separation D « 50 \ ^ , and around station zenith sin 0 ^ 1, so that 
Eq. (33) is 


d«’--5^dCo <34) 

In the present system f^ is 100 kHz, f^ is 100 Hz, and dC the quantization 
of the count is 1. 


'^®ldc„=ieount--2a0'=') radians 


(35) 


Equation (35) expresses the angular error in the liner-of-sight for a one count 
error in the measurement count. Equation (35) expressed in arc seconds is; 

d^lnr . ^--4.1 arc seconds 

‘^’^0 “ 1 count 

and reflects the maximum error in the present Minitrack system due to quanti- 
zation of the measurement. 
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The second source of error considered here is the uncertainty or bias in 
the transmission frequency of the satellite, . (See Appendix A for the insighifi-; 
cant error caused by time varying f^ ). From equation (32) 


- sin 0 d 0 = - 
D 



(37) 


Again baseline separation D ~ 50 \ at zenith K = .0, and at ^ = 85°, K = 4, 


“ *02 radians 

Cq = 999 

■ dft 

^^15 = 85° "0.1 — radians 
Cq = 999 ^ 

K=4 


(38) 


(39) 


d^|g_^o® « 4. 1 (10^) arc seconds 

Cq = 999 t 

. dft 

d(9|-^_85<> » 2.1(10'^) arc seconds 

Cq = 999 t 


(40) 


(41) 


K = 4 

The- restriction of $ in Equations 38-40, is made because the main antenna lobe 
of the fine Minitrack antenna is approximately 10 degrees wide, (See reference 3) . 
Measurement data are taken as the satellite traverses this 10 degree beamwidth. 
Assume that df^ is 1 kHz and that f^ =136 MHz. Then 


^^\\e = 85° “ ■ 

Cg = 999 
A = Ik Hz 


( 42 ) 
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^^l<9 = 9o® ”’~.03 arc sec (43) 

Cq = -999 
A = Ik Hz 

Equations (42) and (43) express the maximum error in the present Minitrack 
system due to an unknown offset of the satellite transmitted frequency (constant 
doppler or bias offset) for the conditions noted. This error source must be 
termed insignificant when viewed against other error sources (also see Appendix 
A for the smaller error caused by doppler rate in f ^.). 

The third and major source of error in the Minitrack system (and any inter- 
ferometer system) are unknown path delays in the measurement system, especially 
the antennas, prior to summing of the signals.- Again refraction .has not been 
made a consideration here. See Figure 2 for an illustration of unknown path 
delay errors. Assume that there exists an unknown delay of r seconds in the 
path from antenna ”A" to the summer. It follows that equation (26) becomes: 


5' = --/-(t - A 

2 77- f 


(44) 


For a constant frequency sinusoid as system input, equation (29) becomes: 


Co 


=_L_r 

2t7 f L 


2 7T D cos ^ 


- 2 77- f . 


(45) 


Define; 


A S = S _ 


(46) 
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At zenith, 0 = 90°; cos ^ = 0 and equation (46) is 


f . = constant 


2 7T f , 


[2w ft t \ 


(47) 


Assume r = 1 nanosecond == 10**® seconds, that is there exists an unknown or 
uncalibrated delay in the antenna ”A" path of 10~9 seconds.4 Then 

seconds 

f ^ = constant = 136 M Hz ^ 

r= 10“® 

From Equation (32) we have 


cos 6 - ^ 
D 



(49) 


- sin 6 dd - — — ^ dS 

D ft 


(50) 


At zenith, and with D = 50 





(51) 


From equation (48) for 


d ~ 90°, = 136 M Hz, T = 10"^ 


^ This is equivalent to a path length change in antenna A, its associated cabling or its electronics 
of 0.3 meters for a velocity of propagtion assumed equal to the speed of light, c. 
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d.5 I a 90° 

f j = constant a 136 M Hz 
T = 10 


A5| 


e - 90 ° 

f ^ = constant = 136 M Hz 


T = 10 


.136 

f 

r 


seconds 


(52) 


and 


f ^ = 136 M Hz = constant 

T-= 10 


.136 

50 


radians 


(53) 


*^^U = 90° =-561 arc seconds 

= 136 M Hz = constant 

T= lO"^ 

This value of 561 arc seconds demonstrates the difficulty in interferometer 
measurement. Flyover calibrations which measure or calibrate the unknown 
path len^h delays in.the system prior to combining are done periodically, an- 
nually or semiannually. To maintain an absolute pointing error of, say 6 arc 
seconds, the path length differences in the antenna,, cables, and electronics prior 
to combining, must remain stable to approximately 10“^^ seconds or an equivalent 
path length difference of 3 millimeters. 


Historically, these Minitrack calibration numbers are stable to the order of 
10 arc seconds. For example, examine graphs 1 and 2 in Reference 5, which 
portray the historical record of the aircraft calibrations at the Winkfield 
Minitrack station from 1961 to the present (20 calibration points determined 
over that period of time). The shifts in the calibration for the antenna, cables, 
etc., for the six month period between August 1969 and February 1969 were as 
follows (4 arc seconds % 1 coxmt) 
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EAST-WEST FINE EQUATORIAL 
NORTH-SOUTH FINE EQUATORIAL 
EAST-WEST FINE POLAR 
NORTH-SOUTH FINE POLAR 


-24 arc seconds 
+ 8 arc seconds 
- 8 arc seconds 
0 arc seconds 


in addition, over the 9 year span of calibration, the follomng maximum spreads 
have been observed for the calibration numbers at the Winkfield station (4 arc 
seconds ~ 1 count): 


EAST-WEST FINE EQUATORIAL 
NORTH-SOUTH FINE EQUATORIAL 
EAST-WEST FINE POLAR 
NORTH-SOUTH FINE POLAR 


96 arc seconds 
60 arc seconds 
52 arc seconds 
48 arc seconds 


These numbers and similar numbers from other Mintrack stations are not to 
be construed as a criticism of the system; they are meant to demonstrate the 
difficulty in establishing pointing accuracies on the order of 10 arc seconds . 
The authors believe that 10 arc seconds is the magnitude of error which should 
be used for system resolution. 


SATELLITE DYNAMICS 

The authors are indebted to Messrs, Simas and Santarpia of GSFC for the 
derivation of the formulas describing the phase of the input signal related to 
satellite dynamics (see Reference 6). ThesC formulas and their pictorial 
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representation are included in the following paragraphs. In Figure 1, Geometry 
of Circular Orbit, a satellite in a circular orbit at a height, h in meters above 
a spherical earth, -passes directly over and parallel to the baseline separation, ‘ 
D in meters, r is the earth radius in meters, e is the Hne-of-sight angle in 
radians with respect to the local horizontal, and e is shown increasing in Figure 
1. The angle a , defined as the angle between the earth radius vector r to the 
station and the vector to the satellite R, is decreasing in Figure 1. The re- 
maining angle, /3, is- decreasing with time. 

a + ;6+ (e+f) = ir (50) 

d -f yS + (9 = 0 (51) 

where the dotted quantities are the time derivatives of the respective angles in 
radians/sec. 

4 = ( 52 ) 

(53) . 

(54) 


.0 + 
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where is the tangential velocity of the satellite along the circular path above 
the earth, is a constant positive quantity. The rotation of the eaai:h has been 
neglected here. From the law of sines 

r _ *~ •*• *^. _ R (55) 

and equation (55) becomes 

r cos,j9 “ R sin jS 


^ = sin~^ cos ^ 

By taking the time derivative of equation (57), one obtains: 

^ = r^sin^ - (58) 

One can verify the correctness of the sign associated with equation (58) in the 
following way; 




and sin0(O<^<7r) 
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are positive quantities. 6 as shown in Figure 1, is always increasing with time, 
so that 6 is always positive. Hence /3 is a negative quantity which physically 
matches the geometrical conditions shown in Figure 1. 

By equating equations (53) and (58), we obtain: 


e = 


i[-(- 


cos d'^ 


(59) 


[‘ - (^“fr 


sin 6 

R 


At ^ = 0 on the horizon, 6 = v^/R, and directly overhead, B -tj /‘I, B = j 

Figure 5 shows 6 , and 0, for the condition of h^ = 100 statute miles, the assumed 
worst case satellite pass. The tangential velocity in meters/second is calculated 
from; 


V 



(60) 


where fj, is the earth’s gravitational parameter, 3.986016 (10^"^ ) m Vsec^ . Time is 
calculated from 


t = 


r + h 



(61) 


with t = 0 seconds at 6 = 0, 



For further calculation, 6 is required. This requires the differentiation with 
respect to time of equation (59). 

r, /r cos /r\2 sin(2<9) • 

(s j ^ 



From the time variation of ^ , it is necessary to calculate a phase variation 
through the system, especially through the postdetection filter shown in Figure 2. 
In Figure 2 and equation (1) we have defined a length 0 where: 

$ = D cos 6 (1) 

and ^ is the difference in path length to the two antennas. Equation (15) ex- 
pressed in terms of D cos 0 is: 

cosja 7T f_. t + 0 (t (t)]| 


20 



which is the postdetection filter output in terms of the changing length D’.cos d » 
It is convenient here to disregard the phase changes introduced by antenna phase 
patterns, but later these effects will need to be considered. Equation (63) with 
the simplification is: 


: |2 77 t + 0 yp Ct)| 


(64) 


Again assume thotip (t) = 2nft; that is, the antenna input is a constant frequency 
sinusoid. In the preceding section, and in Appendix A, it was shown that the sys- 
tem error (error in iine-of-sight angle ^ ) is relatively insensitive to changes 
in f ^ . With this assumption, equation (64) becomes: 

cos [2 77 t 4-^(t)3' 

and expressed in terms of D cos 0 , equation (65) becomes: 




2 77 f D cos & " 


( 66 ) 


Now c - f X and equation (66) is 


; ^2 7T t - cos 


where 


(67) 


^(t) : 


2 77 D cos $ 


( 68 ) 
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and equation (68) states the change in phase of the signal (other than the linear 
change 2 Trf^t) caused by the satellite dynamics. 


From equation (56) 

R 

cos 9 ^ — sin 
r 

and 


(69) 


0(t) = 


2 7T D R . n 
— Sin p 


By taking the time derivative of equation -(70) 


(70) 


<jS) (t) = - ^ ^ — (cos/3)/3 


(71) 


but from equation (53) 


< 5 & (t) = - 5 (cos ^') Uo + 

K ^ V R / 


(72) 


<^(t) 


2 77 D ^.6 cos /S 2 77 D 


(73) 


The time derivative of equation (73) is: 


<!> (t) = 


2 77 D P sin p 


(v^ - B0) 4- ^9cosP 


(74) 


The quantity 0 (t) is also required and a numerical differentiation was performed: 
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(75) 


<f> (t). 


(t = t^) -<^ (t = t^) 


where the time steps -t^ are equal to the time required for 6 to change by 
1/2 degree. This was precise enough for the purpose here. 


Figures (6) and (7) show^, ^Modulo 2 tt ^ D - 50 , and h^ - 

100 statute miles, cb corresponds to equation (19) where the ambiguous 

Modulo 277- 

measurement (to multiples of 2-77) is defined. In Figure 6, 27 t 

6 cycles plotted. In the approximately 37 seconds remaining before zenith, there 
occur 44 more cycles of 2 t 7-. At zenith is changing at a rate (maxi- 

mum) of 15.2 radians per second, or 2.4 cycles per second. From equation (68) 
we have the relationship between^ , the desired result and<^ , the phase function 
of the filter input signal: 


<p (t) = - ^ cos 6 
X. 


( 68 ) 


For D equal to 50 A. ^ 


<f> - - lao TT cos 6 

(76) 

d 4> - 100 ‘TT sin ^ d ^ 

(77) 


Near zenith, sin ^ = 1, and 


d 


d 6 (radians') ^ 


100 TT 


(78) 



where d<^ is expressed in radians. For d (9 in arc seconds: 


d 6^ (arc seconds) » 657 d(^ 

SYSTEM DESIGN CONSIDERATIONS 

Some comments about the signal dynamics and electronic system design are 
needed at this point. These comments are based upon the ’’worst case" assumption 
of a 100 statute mile high satellite passing directly over the station. Refer back 
to Figure 2. In this block diagram much has been eliminated for the sake of 
simplicity. For example, the conversion of the summed antenna A + B signals to 
an intermediate frequency and the filtering (predetection filtering) is not shown. 
No attempt was made to investigate the possible replacement of this passive 
filter and the amplitude detector with a coherent detector (phaselock loop).. 

It is assumed here that for reasons such as ease of operation, the predetection 
filtering will remain a passive filter with a nominal 10 kHz bandwidth and this 
filtering is assumed to have no effect on the signal dynamics. For the post- 
detection filter, the signal dynamics are described by equation (68): 

cos 0 (t) (68) 

and the postdetection filter, whether active or passive must be designed for 
these signal dynamics. The present filter is a passive, 10 Hz (nominal) bandwidth 
filter, which is sufficient for the "worst case "signal dynamics shown in Figures 6 
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and 7 . Time delay through the passive filter, as long as it remains a constant 
over the total bandwidth (not feasible), and over temperature and other variables, 
would cause no error in the system measurement,® 

Here it will be considered that a phaselock loop filter could be built which 
matches as closely as possible the signal dynamics. The next section of this 
report will examine the phase error of a phaselock loop filter for various orders 
and bandwidth, and how the small but finite phase error impacts system error. 

PHASELOCK LOOP RESPONSE TO SATELLITE DYNAMICS 

In Figures 6 and 7, it may be seen that even for the "worst case" satellite 
orbit (100 statute miles directly overhead), the magnitude of the higher order 
derivatives otcf> (t) diminish rapidly with order. This immediately suggests ex- 
pansion of <^(t) in a Taylor series of ^.(t) around any time point t = t^. Thus 
<^(t) can be expressed in the neighborhood of as: 


(t) - 4 ) 4- 0 (t^) [ t , 






tj-' 


* (V 

•3 ! 


ft - tnP 


+ higher order terms. 


(80) 


where the higher order terms are considered negligible here. The value 0(t^), 

^(t^) and <^(tQ) may be read direetiyTrom Figures 6 and 7 at any arbitrary 

time point, t_. For example, at ^ - 85°, and t == 184.4 seconds: 

0 0 


Mr. W. Rice of GSFC is presently investigating the variation of time deloy and its effect upon 
present system error. 
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4>{t^')=-2.2 radians®* 


(81) 


^(0 = 15.1 (82) 

^ sec 

^ (t„) = . 19 (83) 

(sec) 2 

(t„) = -0.10 (84) 

(sec)^ 

For the purpose of solving for a phaselock loop response around the time, t = 
184.4 seconds, the input to the phaselock loop can be represented adequately by: 


4 >.^ (t)« - 2.2 + 15.1 t + 


0 . 


12 t2-^t3 


(85) 


where 0i„(t) is the input to a phaselock loop at, and around t = 184.4 seconds 
for the assumed "worst case" satellite orbit, h^ = 100 statute miles. In equation 
(85), t = 0 at to = 184.4 seconds. 

Table I lists the equations for solution of the error response of an acceptable 
third order loop (with specific characteristics defined by Reference 7) , (t) 

describes the error response 




( 86 ) 


® The system cannot recognize multiples of 2 7T radians, so that <^(to) 2tt 
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TABLE I 


Equations for Error Response-Mallinckrodt Third Order Loop 


I 0 >. 


S +— I (S 


a ('£ 21 ') = 

^Ksec) 0.743 


^in 


(s + j) (S 


is) =\t * ^ t /1 

s*^. 


:T 

J / ^ \ 

|s (s + Oi^) 


(s) =4>__J. ^<b ^ 4 - ^ ^ 

(s+-^)(s+«)2 (s+^)(s+«„)2 (^s + (S i- 0,^ 


(s 4. _j (S X CtJ^)2 
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TABLE I (Continued) 
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exp 


/_%‘\r 32 ^' 

4 /L 3„3 + <„2 

n n 


4 <?b 



-f 


rl4c^ 80 40 

exp ( - o> t) - _ + _ 4- 8 0 N. 


L 0 ?- 


f20 20 ' 1 20 

exp (-co^ t) 4-40-4 cy^0 + (I~7) 




where form of equation (85) with constant coefficients for t, t^, 

and t^, and (t) is the output phase corresponding to the input The func- 

tions of "s’ 'shown in the table are the Laplace transforms of the input function 
and the Laplace transforms of the phaselock loop component walues. 

The solution for the error (t) assumes that the phaselock loop has no , 
initial information about the input (t) , and that the input conditions remain 
over all time of the form: 

(t) =Kj +K 2 t +K 3 t 2 +K 4 t3 (87) 

where K^, K^, K^, and are constants. 

Even with these restrictions, the application of these methods will provide 
insight into the phaselock loop design requirements. Admittedly by choosing 
the Mallinckrodt form of the 3rd order loop, further restriction has been placed 
upon a 3 rd order loop design, but the prior operational experience gained with 
this 3rd order loop (for example, the Goddard Bange and Eange Rate Tracking 
System) and its demonstrated stability, make it an acceptable potential design. 

The one parameter left undefined in Table I is the phaselock loop bandwidth, B^. 
This is the standard symbol for the (one-sided) loop 'noise bandwidth, and for the 
purpose here will be fixed at either 10, 5, or 3 Hz, and- in a few instances at 20 Hz. 

In Table 11, the equations for a second order phaselock loop are given, also 
with as a variable parameter. The damping factor is fixed at 1/V2. The 
equations are given similarly in terms of constant values for <p and <^. 
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TABLE n 



Figure 8 portrays the absolute value of phase error, | 0^ (t) | , as a function 
of time after the initially quiescent third order loop attempts signal acquisition. 

In Figure 8, the satellite is at an elevation angle (interferometer angle) of 85 
degrees at TIME equal to zero. The abscissa value is labelled TIME to differentiate 
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it froiB tile abscissa labelled t in Figures ( 6 ) and (7). At t = 184,4 seconds in 
Figs. 6 and 7, TIME = 0 in Fig. 8 . For TIME greater than zero in Fig. 8 , the 
values of 3 <p»4>» and are assumed to remain constant at the values for 

t = 184.4 seconds in Figs. 6 and 7. Absolute magnitude of ^ was chosen for the 
ordinate of Fig. 8 to avoidpositive and negative values for < 56 ^ . By choosing | |- 

as the dependent variable, the plot, becomes. appro:sim ate ly the envelope of the 
phase error, except at TIME very large, where it approaches 2^/3 . Since 

only a finite number of values of ITME were used in the solution of | <^^ | , (TIME - 
0, 1 / 4 , 1 / 2 , 3/4, 1, 2, 3, 4, 10, 20, 100), the graph has a jagged appearance. This 
is caused by the oscillatory nature of 0^ for TIME less than approximately 2 
seconds. Fig. 8 was constructed by connecting the computed values (0^| at TIME = 
0, 1 / 4 , 1 / 2 , etc. with straight lines. For small values of | 0 |. corresponding to 
large values of TIME, this method of presentation is satisfactory. The reader 
should note the variation of TIME scale along the abscissa. 

The four horizontal lines in Fig. 8 have the following significance: the top 
line shows the approximate limit for computation of error from the linear equa- 
tions of Table I. Above that top line, values of 0^ must be qualified. For example, 
in Fig. 8 , at TIME == 0, a phase exmor 1 0 ^ 1 of 2.2 radians is plotted and this value 
lies above the approximate limit for loop acquisition. By examining tbe linear 
equation (1-7) of Table I, one sees at TIME = zero that 




Mod 2w 


(77) 
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In Fig. 9, 27 t plotted for t between 184.2 and 186.2 seconds. At t = 184.4 

seconds {6 = 85°) is -2.24 radians. From equation (77), at the time of 

attempted loop acquisition of 184.4 seconds, the phase error calculated from the 
linearized phaselock loop equation would be -2.24 radians. Of course, the linear 
equations do not hold for this condition. Instead a series of events occurs: 

(1) An error voltage is generated in the phase detector and passed through 
the loop filter to drive the voltage controlled oscillator (VCO). 

(2) The VCO phase and phase rates will not match the incoming signal since 
the initial error volt^e was generated by the phase detector operating . 
outside its linear range. 

(3) The loop will continue to generate "useless" VCO control information 
until the difference between the input phase and the VCO phase is less 
than approximately one radian,, so that linear operation of the loop is 
begun. From that time onward, the loop will tend to reduce the difference 
between the input phase and the VCO phase (phaselock) so long as the 
dynamics, of this signal are within the tracking capabihly of the loop. 

In Fig. 8, at 0 = 85°, one must assume that the loop is incapable of acquiring 
rapidly. Fortunately from an examination of Fig. 9, a mere 0.1 second later the 
phase input to the loop has decreased to less than 1 radian, so that linear opera- 
tion of the phase detector can begin. These preceding comments are best 
illustrated in Fig. 10. In this figure, acquisition is begun at a time when 
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IS zero \{0 ~ 90°). The phase error, <p^ , at TIME = 0 is very near zero. Except 
for a loop bandwidth of 3 Hz which is incapable of handling the high rates of 
change of phase, the phase error with time remains less than 1 radian (loop 
lock). The explanation for the' 3 Hz bandwidth difficulty can be seen readily from 
Fig. 7. At & - 90°, 0 is equal to 15.2 radians /second or equivalently 2.4 cycles 
per second. Thus the input signal to the phaselock loop is offset in frequency 
by 2.4 Hz, and almost outside an assumed loop noise bandwidth of 3 Hz. 

The other horizontal lines in Figs. 9 and 10 show the error in the inter- 
ferometer angle 6^ in terms of . The values of 0^ .were calculated by means 
of equation (79) 


d 6 (arc seconds) = 0^ ~ 657 dcp = 657 

Now if one assumes that the error contribution from the phaselock loop 
filtering should be less than l/lO of expected system resolution, i.e., less than 
1 arc second (^ , ^ 1 arc second), it can be seen that for a bandwidth of 5 Hz a 
time of approximately 4 seconds is required after acquisition is attempted 
before the interferometer error 9 is less than 1 arc second. This condition 

te 

applies for the satellite passing zenith or near it. Similarly 2 seconds are re- 
quired for the 10 Hz bandwidth and 1 second for a 20 Hz bandwidth. These 
seemingly small lengths of time for signal acquisition^ are significant in terms 

^ Signal acquisition is defined here as that point in time at which 9 ^ is less than or equal to 1 arc 
second. 
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ofthe "worst case" orbit assumptions. For a 100 statute mile satellite height, the 
satellite remains within ±5 degrees of zenith for only 3.6 seconds. The total 
pass, visibility from horizon to horizon, occurs within a span of 6.2 minutes. 

The reader may argue that it is not required to attempt acquisition at or near 
zenith. He may contend that acquisition should occur at the horizon, and by the 
time the satellite crosses station zenith 3.1 minutes later, the phaselock loop 
error should be negligible (less than 1 are second for 6 ^ ) . The authors offer 
the- following refutation. 

Reference 3 describes the amplitude patterns of the Minitrack "fine" antenna. 
The amplitude pattern is seen to be fan shaped and this fan shape serves a useful 
purpose. In general, the satellite passing over the station does not pass exactly 
over the center of the station. For example, it will pass north or south of the 
EAST-WEST FINE EQUATORIAL antenna pair. The fan portion of the beam 
extends approximately 38® to the north and 38® to the south to cover this eventu- 
ality, and the satellite as it makes it traverse from west horizon to east horizon, 
passes through antenna sidelobes, except at approximately ±5 degrees of the 
station meridian. In theory, as the spacecraft signal encounters an antenna 
amplitude null, such as at ±5 degrees off station meridian, it sees a change of 
phase of 180° caused entirely by the transition from the first antenna sidelobe 
to the main antenna lobe. It is useful here to refer back to equation (16) 

= {i// (t - At) - i/; (t) + V (t - At)] -U (t)l } (16) 
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U and V have been described as the phase changes caused by antemias A and B 
in Fig. 2, If the two antenna phase patterns were identical and if At were equal 
to zero, V would be equal to U for any satellite dynamics and there would be no 
contribution to <^{t) from antenna "A" and "B" phase patterns. In practice, the 
phase patterns are unequal and A t is only zero at station zenith (or station 
meridian for our example) . For the phaselock loop which has acquired at the 
horizon, as it observes the- signal passing from antenna lobe to antenna lobe, 
it will see a time varying phase perturbation caused by the Minitrack antenna 
pair. The authors contend' that signal acquisition at the horizon is no better than 
signal acquisition at ^ = 85 degrees, since the phaselock loop will require 
time to reduce 6^ to less than 1 arc second for each phase transient it sees at 
its input, whether it is caused by signal dynamics, or antenna phase patterns. s 

Figs. (8) and (10) covered the conditions for a Mallinckrodt third order loop 
•serving as the^, postdetection filter. The next figures will show operation with 
a second order loop (damping factor = 1/V2) for the same conditions as for 
Figs. 8 and 10. 

Figure 11 has acquisition attempted at 0 = 85*^. The comments for exceeding 
linear operation of the loop also apply here. There will occur some small but 
finite time after which the loop will assume linear operation. For the second 
order loop, and a non-zero 4 > , it may be observed from Equation (II- 8) in Table II 

® No quantitative data exists to verify this contention. A source of data could be a Minitrack station 
calibration wherein the calibration'plane flies a pass parallel to the baseline rather than normal 
to it. 
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that will increase linearly with time. At a very large.time, will be very 
large. Recall, however, the assumptions made about the input to the loop. by 
assumption is to remain constant at all time for the,value encountered at ^ = 

'85°. But for the "worst case" orbit pictured in Fig. 7, ^ is significant only over 
±20 seconds about station zenith. For the 10 Hz loop bandwidth of Fig. 11, it 
requires 100 seconds of. constant </> to reach an interferometer angle error, ^ 

e 

of 3 arc seconds, but physically cannot remain a constant non-zero value over 
that length of time. 

Acquisition times for the 10 and 5 Hz bandwidth are on the order of l/2 and 
1 second. In Fig. 12, acquisition attempted, at 0 = 90°, ^tt zero, and the 
acquisition time for the 10 and 5 Hz bandwidths remain near l/2 and 1 second. 

The 3 Hz bandwidth is obviously too narrow as evidenced by its transient error 
exceeding linear loop operation. Phase error in Fig. 12 also is seen to increase 
linearly with time, but the same qualifications expressed in the preceding para- 
graph apply here as well. Third derivative of phase cannot remain significant 
over 100 seconds of time. 

Figures 13 and 14 illustrate third order loop response for a 500 statute mile 
circular orbit. For this orbit the satellite remains within ±5 degrees of station 
zenith for 19 seconds. The entire pass from horizon to horizon occurs within a 
time span of 15.4 minutes. Acquisition times have not decreased notably over 
the 100 statute mile "worst case" orbit (Figs. 8 and 10). Here the 10 Hz loop 
has an acquisition time of approximately 1.5 seconds, the 5 Hz loop approximately 
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3 seconds, and the 3 Hz loop approximately 4.3 seconds. In Fig. 13 the values 
of above 1 radian must be qualified as in the previous figures. Note in Fig. 14 
that a 3 Hz bandwidth does not exceed linear loop operation. For the 500 statute 
mile orbit at 0 =90°, is a maximum and is equal to 2.9 radians/second, or 
equivalently a frequency offset of 0.5 Hz. The 3 Hz bandwidth, for these con- 
ditions, does not exceed a phase error of 0.22 radians, but does have a long ac- 
quisition time, 4.3 seconds. 

Figures 15 and 16 are for the condition of second order loop and the 500 
statute mile satellite height, at ^ = 85 and ^ = 90 degrees. The satellite dy- 
namics are identical for Figs. 13, 14, 15 and 16. Acquisition times are approxi- 
mately 1/2, 1, and 2 seconds for the 10, 5, and 3 Hz bandwidths. Since a constant 
value of 4> is assumed at time of attempted acquisition, increases with time. 
For example, in Fig. 15, the included table gives values of for large values 
of TIME; at TIME = 10,000 seconds, and = 10 Hz, cj>^ = 2,1 (10"^) radians 
and - 1 arc second. It must be repeated here that a constant <f) over 10,000 
seconds is not physically realizable. The inclusion of a constant <f> was made 
to show the unimportance of this source of error for the second order loop. 

RECOMMENDATION FOR LOOP BANDWIDTH AND ORDER 

The preceding material has defined sufficient constraints in order that a 
suitable design may be made for the postdetection filter of the Minitrack system. 
The assumptions made in this definition are repeated here: 
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• The "worst case" conditions occur for a circular orbit of 100 statute mile 
height passing directly over the Minitrack: station. In this situation, satellite 
dynamics define maximum rates of 4 > , the input to the postdetection filter. 

o No design changes are anticipated in the predetection filter and the 
amplitude detector, to reduce or increase the rates of into the postdetection 
filter. 

• Phase transients caused by input signal entrance through antenna side- 
lobes and nulls prevent "accurate" data taking near station horizon,® and force 
acquisition and "accurate" data taking at ±5 degrees about station meridian 
assuming an east- west pass. 

o An error contribution of 1 arc second from the postdetection filtering 
is desirable and practical. This 1 arc second is l/lO of expected overall system 
resolution (10 arc seconds). 

• Acquisition time (reduction of $ to 1 arc second) can be calculated from 
linear phaselock theory and Taylor series expansion of <^(t), the*phase function 
into the phaselock loop. 

• 4> cannot remain an important Input to the phase lock loop over long 
periods of time. See Fig. 7 for illustration. 

With all the preceding qualifications listed, it is recommended that a 10 Hz 

bandwidth, second order phaselock loop be utilized as the postdetection filter for the 

9 Calibration of antenna pairs is done no lower than ±50 degrees zenith, so that data taken outside 
this calibration zone is questionable (see Ref. 1). 
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Minitrack system. Acquisition times for the ’*worst case" orbit conditions are on 
the order of 1/2 second, and the error remains well below 1 arc second in pointing 
angle for any reasonably assumed ^ . A 3 Hz bandwidth second, order loop is 
inadequate for following the large ^ of the "worst case" orbit, and cannot be con- 
sidered acceptable. A 5 Hz, second order loop doubles the acquisition time to 
approximately 1 second, or approximately 25% of the total visibility times be- 
tween ±5 degrees of station zenith for the 100 statute mile satellite height. 

The acquisition times for the third order Mallinckrodt loop are large frac- 
tions of the visibility times between ±5 degrees of station zenith. The advantage 
of the 3rd order loop lies in the reduction of error caused by ^ , but <for the "worst 
case" orbit assumedhere, 0 is not an important factor in system error. For this 
reason, the 3rd order loop described in Table I is not recommended for 
implementation. 

The reader will notice' that no mention of signal to noise ratio has been made 
to this point. Except for the assumption of a total system resolution of 10 arc 
seconds, based upon repeatability of aircraft calibrations, no error contribution 
from noise has been calculated. System design must be based upon anticipated 
signal dynamics, and when the design is completed, parameters in the link calcu- 
lation which affect noise errors in the system must be sized accordingly. A 
sample calculation will be given here in explanation of this important point. 
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The .ppstdetection filter is here required to be a 10 Hz, second order loop 
filter from the expected satellite dynamics. Preceding the postdetection filter is 
an amplitude detector, with a 10 kHz filter preceding the amplitude detector. 

System noise figure is assumed to be 4 db, and mtenna gain at ±5 degrees about 
zenith is approximately 15 db. If the pointing error from signal to noise consider a- 
tipn is required to be 1 arc second, then from Equation (88) 

= 1 . 5 (10“^) radians 
From linear phaselock loop theory 

. <#., = 1 . 5 •( 10 - 3 ), 

where S/N is the signal to noise power, ratio into the loop, and 2 4> is the 

noise power, with 0 as noise power density. Since our recommended = 10 Hz: 

2. 22 (lOS) =^(53.5 db) (91) 

4-44 (10®)=^ (66. 5 db^^.Hz) 

where S/4> is the required signal to noise power density ratio for a pointing error 


of 1 arc second. 



Table m gives sample calculations for required signal input for .a pointing 
error of 1 and 10 arc seconds. Antenna gain is given for ±5 degrees about zenith. 

At antenna nulls, a reduction of signal from 20 to 30 db can be expected. The ampli- 
tude detector then will see signal to noise ratios of less than one for the inputs* of 
Table m, and signal to noise ratio at the detector output will be equal to the 
square of the input signal to noise ratio. Consider here that data is taken ; 
only in the main lobe, i.e. antenna gain = +15 db. For a* 1,000 statute mile range 
to the satelhte at a frequency of 136 MHz, the space loss is 139.5 db. It follows 
that atthe satellite an effective radiated power (ERP) of +21.0 dbm is required 
for a 1 arc second noise error, and +1.0 dbm is required for a 10 arc second 
noise error. This ERP is for carrier signal only, if modulation occurs during 
the Minitrack measurement, the total satellite ERP must be increased to account 
for the power in the modulation sidebands. 

Now assume that for the signal inputs of Table III, the satellite is located 
at an antenna null of -30 db. Table IV lists the calculations for this situation, 
and it is apparent that drastic changes' in pointing error occur while passing 
through a 30 db antenna null. For the input signal level corresponding to the 1 
arc second erroi* while in*the main loDe oi the antenna, an error of 46.4 arc 
seconds occurs at the assumed 30 db antenna null. For the input signal level 
corresponding to the 10 arc second error, no valid measurement' is obtained 
during the antenna null (probable phaselock loop unlock); 
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TABLE m 


Signal to Noise Calculations - Satellite in Main Antenna Lobe 


6=1 Arc Second 

e 

Required Signal Level Into Antenna 

-118.5 dbm 

Minitrack Antenna Gain (Main lobe only) 

+15 db 

Required Signal Level Into System 

-103.5 dbm 

Noise Power in 10 kHz Predetection 


Bandwidth (System NF = 4 dB) 

-130 dbm 

S/N Ratio Before Amplitude Detector 

+26.5 db 

S/N Ratio After Postdetection Filter (2Bj^ =20 Hz) 

+53.5 db 

Pointing Error 

1 arc second 

e =10 Arc Second 

e 

Required Single Level Into Antenna 
Minitrack Antenna Gain (Main lobe only) 

-138.5 dl^ni 
+15 db '-' 

Required Signal Level Into System 

-123.5 dbm 

Noise Power in 10 kHz Predetection 


Bandwidth (System NF = 4 dB) 

-130 dbm 

S/N Ratio Before Amplitude Detector^® 

+6.5 db 

S/N Ratio After 'Postdetection Filter (2Bj^ = 20 Hz) 

+33.5 db 

Pointing error 

10 arc second 


With signal to noise ratios of greater than one into the amplitude detector, output signal to 
noise ratio is equal to input signal to noise ratio. 
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.'TABLE' IV 


Signal to Noise Calculations - Satellite in Antenna Null 
Signal Input Levels Same as Table in 


Assumed Signal Level into Antenna (From Table III) 

-118.5 dbm 

Minitrack Antenna Gain (Antenna Null) 

-15 db 

Signal Level into System 

-133.5 dbm 

Noise Power in 10 kHz Predetection Bandwidth 


(System NF = 4 db) 

-130 dbm 

S/N Ratio Before Amplitude Detector^ ^ 

-3.5 db 

S/N Ratio After Amplitude Detector (10 kHz BW) 

-7.0 db 

Pointing error (antenna null) 

46.4 arc seconds 

Pointing error (antenna main lobe) 

1 arc second 

Assumed Signal Level into Antenna (From Table III) 

-138.5 dbm 

Minitrack Antenna Gain (Antenna Null) 

-15 db 

Signal Level into System 

-153.5 dbm 

Noise Power in 10 kHz Predetection Bandwidth 


(System NF = 4 db) 

-130 dbm 

S/N Ratio Before Amplitude Detector 

-23.5 db 

S/N Ratio After Amplitude Detector (10 kHz BW)1 

Below Loop 

S/N Ratio After Postdetection (2B^^ = 20 Hz) J 

Threshold 

Pointing error (Antenna Null) 

no measurement 

Pointing error (Antenna Main Lobe) 

10 arc seconds 


With signal to noise ratios of less than one into the amplitude detector, output signal to noise 
ratio is equal to the square of the input signal to noise ratio. 

44 




In order to assure at least a, 10 .arc second error during the. 30 db antenna 
null, the required signal level input to the antenna would be -108.5 dbm and the 
equivalent satellite ERP required at a range of 1000 statute miles would be +31.0 
dbm or slightly over 1 watt. 

The repetition here is for emphasis. Signal dynamics define system param- 
eters such as postdetection filter maximum bandwidth and order of loop. Doppler 
shift and phase slope characteristics define predetection filter characteristics. 
Signal to noise constraints are set by maximum satellite ERP. System resolu- 
tion is set independent of all the foregoing, and is strictly a function of time 
delay calibration and time delay variation. 

Still one additional parameter remains to be chosen, and that is system data 
rate. The postdetection filter has been described here as a 10 Hz, second order 
loop, and there exists a reasonable data rate commensurate with this filter type. 
The following sections will concern this optimum data rate selection. 

THE "GAIN” IN SIGNAL TO NOISE RATIO DUE TO LEAST SQUARES 

SMOOTHING OF MINITRACK DATA AT THE OUTPUT OF THE 

POSTDETECTION FILTER 

As previously mentioned, the orbital dynamics of the satellite determine 
the specific type of design for the postdetection filter. Once the design has been 
fixed, the method or manner in which the data is processed at the output of the 
filter determines the sampling rate. 
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In this section, we will investigate the effect of least squares smoothing of 
Minitrack data which has been correlated by the postdetection filter. It is 
assumed that data with white noise having zero mean and luiit standard deviation 
is input to the filter. At the output, the data is smoothed with a least squares 
straight line (it is shown later that a straight line provides a sufficiently good 
fit to the data for purposes of this analysis) and the value of the straight line at 
the midpoint is used as the "value for the data" at that particujar instant of time. 
The standard deviation of the straight line at the midpoint is determined and the 
ratio of the standard deviation of the input noise to the standard deviation of the 
straight line at the midpoint expressed in decibels is used as a measure of the . 
improvement or "gain" in signal to noise ratio due to the least squares smoothing 
process. Three types of filters are investigated — single pole, double pole, and 
ideal (infinitely many poles) with 3-db bandwidths of 3Hz, 5Hz, and lOHz. Inas- 
much as the transfer function of a phaselock loop can be described by an equiva- 
lent electronic circuit, (references), the results of this I analysis for a double 
pole filter can be applied directly to the equivalent second order phaselock loop 
recommended above (the present Minitrack system has a double pole filter with a 
nominal bandwidth of lOHz) . 

In this analysis, we will assume the same "worst case" orbital dynamics as 
assumed in the previous sections hold; that is, the satellite is in a 100 statute 
mile circular orbit, the pass is directly overhead, and the visibility arc is 10° 
(±5° about zenith). In addition it is assumed that the satellite is' acquired 
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immediately 'after entering the beamwidth of the Minitrack' antenna; Although 
there is some finite -delay due to acquisition, this delay is shown to be- negligible 
in its effect upon the "gain.” (see Figure 29.) 

From the orbital dynamics 'and geometry indicated in Figures 1 and 17, the 
"total visible time" that the satellite is in the beamwidth of the antenna can 
be calculated 

T =1^ (93) 

s 

where is the angular rate of the satellite, and is equal to a . 


• CO =1/ 

^ Kt 


+ h 


in radians/second and a is the central angle in radians 


( 9 ^)' 


in-{: 


COS ^ r sin 6 + Vr^ sin^ 6* +. h'^ + 2 r h ] 


r -f h 


(95) 


and 

/Lt= earth's gravitational constant 
r = earth’s radius 

= satellite height above the Minitrack station 
d = interferometer angle. 
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From the above, for a 10° overhead arc and a satellite in a 100 statute mile 
orbit and a 500 statute mile orbit T* =3.6 seconds and =18.9 seconds 
respectively. 

The total number of points N during time.T^ is 



r 


where r = time between samples 

In general, the flow of information is as shown in Figure 18. 

For purposes of this analysis a least squares straight line evaluated at the 
center point provides a good fit to the data. This can be seen in the following sense. 
From Fig. 1, the measurement of phase delay is proportional to cos 6 where d 
is the interferometer angle. Figure 19 shows a least squares straight line fit 
to sin (90 - ^ ) from +5° to -5° where it is compared to values of sin (90 
It can be seen that both curves pass through zero. Therefore in the absence of 
any systematic and random errors the data and the fit to the data will agree at 
the midpoint. It is in this sense that the straight line fit is sufficient for the 
purpose of this analysis. Furthermore, it is felt that the same conclusions can. 
be drawn with a higher ordered polynomial fitted to the data over the same arc 
lei^h. 

In Fig. 18, assume that white noise (uncorrelated) is input to the filter and 
that (2n + 1) data points are available for processing (equally spaced in time 
from - n to +n). At the filter output the (2n + 1) data points (by using 2n + 1 points 
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symmetrically' spaced' about zero'the mathematics Is-ihade simpler); are fit' with 
a least squares' straight line and the' value of the straight line at the center ^point 
e^r'adted as the value for the data at that time.' * The question then- arises,' ‘"-What 
*gaih*ds achieved at the filter output due to least squares smoothing?"' 

It will be seen that the "gain" achieved is a function of the number of data 
points, the type of filter (single pole, double pole, ideal), the bandwidth of the 
filter and the sampling rate (l/r).' The number of points durihg a total time T^ 
is given by equation (96)-. For' computation purposes, the largest number of 
odd points in time T^ was taken for N = 2n h- 1. 


The autocorrelation functions for the single pole, double pole, and ide^ filters 
are respectively (references 9 and 10) 

Single Pole 

P, (97) 

‘ ' ■ i 

where = 3 db angular cutoff frequency of the filter. 

Double Pole 


o> h w h\ oj h 

• C . C \ ^5 C • 

1 + COS 1 e ^ ^ 

fl ! fl 

and 



Ideal Filter. 


sin (i h) 

(i,% h) 
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-.(•98) 


Pi. = 


(99) 



The autocorrelation functions above give the correlation between po'ints - spaced 
0, h, 2h, . . . , (2n h) units apart. It can be seen that by decreasing the bandwidth 
co^ of the filter, the correlation is increased. In equation (96) r is the sampling 
interval which has the same meaning as h in equation (97), (98), and (99). Thus, 
by increasing the sampling rate (or equivalently decreasing the sampling interval) 
the correlation is also increased. 


From the theory of least squares (Appendix B and reference 11) the uncertainty 
of a straight line at the ith data point is given by (matrix notation) 

cr_? = GA^ A)"*^ A^ Q A (A^ A)”^ u (100) 

where 

= (1 i £2 . . . . £2") 

A is the matrix of normal equations, i.e.. 



and Q is the covariance matrix of the input noise (at the filter oulputs). For 
a single pole filter ' 
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where cr^ is the staadard deviation of the white noise at the filter input. 

For a single pole filter, the computations are relatively simple and the 
mean square error at the midpoint of the straight line (i = 0) is given by 



where cr - i. 
e 

For the doiible pole filter and ideal filter, the' computation for the determina- 
tion of orl is rather involved, since it involves the summation of p, and p^ 
yo . 

from 1 to (2n + 1). Therefore, these calculations were done by computer. 

Results of all the computations can be seen in Figs. 20 through 32 where the 
"gain" due to least squares filtering (smoothing) is expressed in decibels, i.e. 


Gain = 20 log, (103) 

where cr = i. 

In each figure, the gain is plotted versus the sampling interval (in seconds). 
Figure 20 shows the effect of smoothing of uncorrelated noise for a 100 statute 
mile overhead orbit (total "visible" time = 3.6 seconds ±5° of zenith). 
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Figures 21 through 26 are for a single pole filter, tlie 'firs't three curves for 
a 100 statute mile orbit and the following three for a 500 statute mile orbit 
(3 Hz, 5 Hz, and 10 Hz Bandwidth). The next three figures (27, 28, and 29) 
pertain to a double pole filter, 100 statute mile orbit and bandwidths of 3 Hz, 

5 Hz, and 10 Hz. Figures 30 through 32 are for an ideal filter, 100 statute mile 
orbit, and bandwidths of 3 Hz, 5 Hz, and 10 Hz. 

It should be noted that the curves are jagged due to the fact that the same odd 
number of points are used for a range of sampling intervals. 

From inspection of Figure 29, it can be seen that the effect of a loss of 0.4 
seconds due to acquisition time is negligible as far as its effect upon the "gain" 
due to smoothing. 

DISCUSSION AND RECOMMENDATION FOR SAMPLING RATE 

From the curves in Figs. 20 through 32, it can be seen that the effect of 
correlation due to the postdetection filter is to reduce the gain achieved by the 
least squares smoothing process. For example, in Fig. 20, the gain due to least 
squares smoothing of white (uncorrelated) noise for a sampling interval of .01 
seconds and a 100 statute mile circular orbit (10° zenith arc) is 25.6 db. Under 
the same orbital dynamics conditions and sampling rate for a 10 Hz single pole, 

10 Hz double pole and 10 Hz ideal filter (infinitely many poles) the gains are 
respectively 20.4 db, 19.2 db, and 18.6 db (Figures 21 j 27, and 30). 

It can also be seen that the higher the sampling interval (the lower the 
sampling rate) the less the correlation between sample points and therefore the 
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effectof'the filter* in this'case approaches that of '-white ridise/ Inspection -of = 
Figs. 21, 22,' 23; and 27 through 32, shows -that the i'gain-** is on the-ofder of' 5 db 
for sampling intervals from 0,73 seconds to’-l seconds (compare: with Fig. -20) 

It is important to mention again that the choice of a particular type of post 
detection filter (e.g. single pole, double pole, ideal) as well as the filter band- 
width should not be made on the basis of least squares smoothing, but should 
depend upon requirements imposed on the system by- satellite orbital dynamics 
as -discussed previously. For example, referring to Figs.- 23, 29, and 32:{100' 
statute mile orbit, 10 Hz single pole, double pole and ideal filter) for a sampling 
interval of 0.01 seconds (100 samples/second rate) the gains are respectively* 
20.4 db, 19.2 db, and 18.6 db; for a 0.05 second interval the gains are 18.2 db, 

17.7 db, and 18.5 db; while for 0.1 seconds the respective gains are 15.4 db, 

15.6 db, and 15.4 db. Thus for a fixed bandwidth and sample rate, much higher 
than the bandwidth, there is little difference among the 3 filter designs -less than 
a db for sampling intervals in the range of 0.05 to 0.1 seconds. 

A question which arises is "For a' fixed filter design such as a double pole* 
filter, can anything be gained’ by increasing the bandwidth? 

At first glance, it would seem that this is the case. However, whenever the 
bandwidth is increased, the noise power increases in direct proportion (reference 
9). This is shown* in Fig. 33, which is -a plot of the degradation in db due to an 
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increased bandwidth. As an example, doubling the bandwidth results in a loss 
of 3 db. Referring, to Figs. 31 and 32 for a sampling interval from 0.05 to 0.1 
seconds, a gain of approximately 3 db is shown from the curves. Therefore, 
the effective enhancement of the signal to noise ratio due to least squares smooth- 
ing in going from 5 to 10 Hz has been offset by a degradation in noise due to the 
increased bandwidth. 

From the above, it can be concluded that the value of the least squares 
smoothing lies in choosing a range of sampling rates to achieve maximum gain 
for a fixed filter design. Thus for the 10 Hz double pole filter design recommended 
in the previous section, it appears that a sample rate of from 10 to 20/sec will 
provide sufficient ’’gain” while keeping the data handling at a minimum. 

CONCLUSIONS 

A study of the Minitrack System, especially in the area of proposed changes 
to the postdetection filter, was undertaken in support of a potential "renovation" 
of the NASA "workhorse". Proposed here as a suitable replacement for the 
present, double pole, 10 Hz bandwidth, postdetection filter is a 10 Hz noise bandwidth, 
second order phaselock loop filter. This proposed filter adequately meets the 
requirement for rapid data taking under the worst case condition of a satellite 
in a circular orbit 100 statute miles above the Minitrack Antenna pairs. 

Minitrack data rate selection is based upon the postdetection filter design 
and that rate is determined from the following considerations: 



(1) For a sampling rate greater than the bandwidth of the postdetection filter 
the effect of correlated noise (as opposed to white noise) in the Minitrack 
data at the filter output is a reduction in "gain" in signal to noise ratio 
achieved by least squares smoothing. 

(2) The choice of a filter design cannot be made on the basis of the "gain" 
due to smoothing at the filter output for the following reasons: 

(a) For a given bandwidth and sampling rate greater than the bandwidth 
there is little difference in the "gain" achieved by smoothing for 

a single pole, double pole-, and ideal filter (an infinite number of poles). 

(b) For a filter with a fixed nximber of poles and sampling rate much 
greater than the bandwidth although there is an "apparent gain" 
from the smoothing process by increasing the bandwidth, this gain 
is offset by a loss due to the fact that the system noise has in- 
creased proportionate to the increase of bandwidth. 

(3) For the 10 Hz equivalent double pole filter mentioned above, a sampling 
rate in the range of 10 to 20 samples per second provides a sxifficiently 
high "gain" when its output is smoothed while still keeping the data 
handling problem at a minimum level. 
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APPENDIX A 


Negligibility of Spacecraft Frequency Rate in Pointing Error Analysis 

The text of the report considered the error from uncertainty in the space- 
craft transmitted frequency, and for a constant tincertainty in the transmitted 
frequency (equivalent to a constant doppler) of 1 kilohertz, the error was shown 
to be negligible. Here we will show the insignificant contribution to system error 
from expected rate of change of received spacecraft frequency. 

In the assumed ’’worst case” circular orbit, 100 statute mile hei^t, the 
range to the satellite can be calculated from: 


RANGE = Vr 2 + r2 - 2 R r cos a 
Rate of change of range can be calculated as 

RANGE RATE = - ^ sin a 
RANGE 

and the time derivative of range rate can be computed from: 


(A-1) 


(A-2) 


RANGE RATE RATE = — — 
(RANGE) 2 


(RANGE) v^ 


RANGE RATE_ 


(A-3) 
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The corresponding doppler frequency and doppler frequency rate can be calculated 
approximately from: 

DOPPLER FREQUENCY = - (RANGE RATE) ^ (A-4) 

^ c t 

DOPPLER FREQUENCY RATE = K = - CHANGE RATE RATE) ^ 

^ c ^ 

For the worst case orbit, varies from 3.5 kilohertz at station horizon to 300 
Hz at $ equal 85® and 0 at zenith, varies from zero at the horizon to -167 
Hz/sec at $ equal to 85® and -169 Hz/sec at zenith. 

Recall the expression for the measurement 8 where 

,.Co. (A-6) 

fr 2 7T f. 

All symbols are defined in the text of the report. Consider that 

(t) = 2 77 (t) t (A-7) 

that is the received frequency at antenna "A" and "B" is a function of time. 
Express f^(t) as a Taylor series expansion about t = 0, the arbitrary reference 
time at antenna "A": 

^ (t) = (0).+ (0) t + *f ^ (0) t2 + . . . (A-8) 
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The reference time is arbitrary because the system concerns itself strictly with 
the frequency at antenna "A” at one instant of time and the frequency received' 
only At seconds before at antenna "B”. Thus 


(t) =2 7 t (0) t + 2 77- (0) 2 ?r y (0) -f • * 


(A-9) 


Similarly 

■/- (t - At) = 2!T (t - At) (0) + f (0) -ft - At]2 + L^'lt - At]2 + . . .|(A-10) 
and 

[./-(t) -0(t - =|2,7 (0) At + 2,77 (0)At [2,t -At].+ • 

At i{t - At)2 + t2 + t (t‘- At)] . . .+1 

Mod 27T 

We will consider here that f ^ (0) and all higher order derivation of f^ at the ’ 
reference time of zero are zero and equation (A-ll) becomes: 

t^<^t) -^(t AAt)]„^^^^ = 2 7r (0)At+-f^ (0)'A t (2 t - A t)];,^^ , (A-12) 

Mod 1 means the fractional part of the bracket left after division by 1. For example 

(A-13) 
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The first term of the bracket contains the desired Information 


f -(O)ifSif • 


(A-14) 


while the second term carries the error, due to a non-zero f ^ (0) 


K ( 0 ^ 


D cos 6 




(A-15) 


The maximum value of (A-15) occurs for the greatest time difference between 
arrival at the two antennas, where t = At, and equation (A-15) becomes 


K w( 


D cos d\^ 


(A-16) 


The maximum value of f ^ (0) is -169 Hz/second at zenith for the assumed "worst 
case" orbit so that (A-16) is 


The ratio of this undesired term (A-17) to the desired term (A-14) is: 
169 

\ c _ 169 CD cos 6) 

f (0) ^ 


(A-17) 


(A-18) 
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For D = 50 X ^ , A-1'8 becomes 


(169) (50) cos e 

ffTO) ^ 


< 4.5 (10“i3) 


(A-19) 


The pointing error from this type of error source is always less than 2 (10 5) 
arc seconds. 
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APPENDIX B 


The Uncertainty of a Least Squares Straight Line at the Center Point 
For Data Correlated by a Single Pole Filter 

Assume (2n + 1) observed values y , at equally spaced time intervals (nor- 
malized to length 1 without loss of generality) . A first degree polynomial 




(B-iy 


' (i = _ n, - (n - i), , (n - 1), n) 


is then fit to the data by the method of least squares. Equation (B-1) can be v 
written in matrix notation 


Y = AX 


(B-2) 


where 


/'-A /' 


ho\ 


Y = 


. A = 


, X = 


V./ 






From the least squares theory the solution X for X is given by 


(*) 


X (A'^ A)"i A^ Y 


(B-3) 


* A is the transpose of the matrix A and (A^- A)**^ is the inverse of (A^ A), 
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which has the covariance matrix 

P = (A*’ A)~i A"^ Q A (AT A)-i 

where Q is the cov^iance matrix of the noise, on the data. 

The variance of the straight line in (B-1) can be written as 

ct2 =0-2 +2X0- + 0-2 

or in matrix notation 

crl = uT (AT A)-i AT Q A (AT A)-i u 
y i 

where 



(B-4) 


(B-5) 


(B-6) 


and a2 is the variance of a. , o- , is the covariance between and a^ , and 
®o ” 

0-2 is the variance of . 


The autocorrelation function for a single pole filter is from references 8 and 
9 given by 


R(ry = 


“ cu 

e 


(B-7) 


■ where co is the 3-db angular cutoff frequency of the filter. If h is the sampling 

C - ' ' > ' ' y ' > - . ‘ 
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interval then 



- i <w h 


(B-8) 


(i = - n, . . . , n) 

gives the coefficient of correlation -between data points spaced ih saihpling’ inter- 
vals apart. 

For white (uncorrelated) noise in the data at the input to the filter having 
zero mean and unit variance, the covariance matrix Q , at .the output of the filter 
will be 



Performing the necessary matrix qperat'ions indicated in (B-6) and letting 
i — 0 (midpoint) • 




2'e 


12 n - . 


2 e 


(1 


n . 




which gives mean square error of the straight line at the center point. The un- 
certainty of the- straight line -at the center is the square- root of (B-1‘0)' 




66 






DATA MEASURE 
(!N COUNTS)' 
TO PUNCHED TAPE 
RECORDER 

0< Co< 999 


NASA-GSfC-T&DS 

MISSION & TRAJECTORY ANALYSIS DIVISION 

BRANCH SSI DATE 7-14-70 

B Y TV PLOT NO. /2<?f 


Figure 2. Simplified Block Diagram of Minitrack System 
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POSTDETECTION FILTER OUTPUT ' 
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Figure 3. Waveforms In Mlnitrack Detection System 
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Figure 4. Time Relationships of Start/Stop R.ulses For Constant Frequency Input 
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Figure 8. Absolute Value of Phose:Error a^ Loop. Output for Time After 
Acquisition Attempted (^ = 85°; t =184.4 Sec into Pass) 
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Figure 11. Absolute Value of Phase Error at Loop Output for Time After 
Acquisition Attempted {6 =85°/ t = 184.4 Sec into Pass) 
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Figure 13- Absolute Value of Phase error ot Loo’p Output for Ts^e Alter 
Acquis I Hon Attempted {$ “85*, f = 451.3 Sec Into Pass) 
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Figure 15. Absolute Value of Phase Error qt Loop Output for Time After 
Acquisition Attempted (^ = 85°, t =451.3 Sec into Pass) 
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Figure 17. “Total Visibility Arc** For An Overhead Pass 
Note: Not Drawn to Scale 
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Figure 18. Flow of Information 





SIN X AND LEAST SQUARES STRAIGHT LINE 
FIT TO SIN X 



Figure 19. Comparison of Least Squares Straight Line Fit to Sin X with Sin X Evaluated 
At Several Points (X =90 ^ = Interferometer Angle) 
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Figure 20. “Gain” Due to Least Squares Smoothing. Uncorrelated Noise 
100 Statute Mile Orbit (T =3.6 Seconds) 
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Figure 21. “Gain” Due to Least Squares Smoothing. Single Pole Filter, 3 Hz Bandwidth 
100 Statute Mile Orbit (T =3.6 Seconds) 
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Figure 22. "Gain” Due to Least Squares Smoothing. Single Pole Filter 5Hz Bandwidth 
100 Statute Mile Orbit (T = 3.6 Seconds) 
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Figure 23. “Gain” Due to Least Squares Smoothing. Single Pole Filter^ 10 Hz Bandwidth 
■ 100 Statute Mile Orbit (T =’3.6 Seconds) 
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Figure 24. “Gain” Due to Least' Squares Smoothing’. Single' Pole Filter 3 Hz. Bandwidth 
500 Statute Mile Orbit (T= 18.9 Seconds) 
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Figure 25. “Gain'* Due to. Least Squares Smoothing. Single Pole Filter 5Hz Bandwidth 
500 Statute Mile Orbit (T= 18.9 Seconds) 
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Figure 26 h “Gain'* Due to Least Squares Smooth irtgn Single Pole Filter 10 Hz Bandwidth 
500 Statute Wile Orbit (T ™ 18-9 Seconds) 
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Figure 27. “Gain” Due to Least' Squares Smoothing Double Pole Filter 3Hz Bandwidth 
100 Statute Mile Orbit (T = 3.6 Seconds) 
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Figure 29. *‘Gain* Due to Least Squares Smoothing. Double Pole Filter 10 Hz Bandwidth 
100 Statute Mile Orbit (T = 3.6 Seconds) 
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Figure 30. “Gain” Due to Least Squares Smoothing. Ideal Filter 3Hz Bandwidth 
100 Statute Mile Orbit (T = 3.6 Seconds) 
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Figure 31. “Gain" Due to Least Squares Smoothing. Ideal Filter 5 Hz Bandwidth 
100 Statute Mile Orbit (T =3.6 Seconds) 
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GAIN" IN SIGNAL TO NOISE RATIO (IN DECIBELS) 
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Figure 32. “Gain’* Due to Least Squares Smoothing. Ideal Filter lOHz Bandwidth 
100 Statute Mile Orbit (T =3.6 Seconds) 
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Figure 33. Noise Degradation. in DB Due to Increased Bandwidth 




